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Abstract 



o 
O 

^ ■ We study the impact of a running index at on the spectrum of relic gravitational waves 

^ . (RGWs) over the whole range of frequency (10^^^ ~ 10^°) Hz and reveal its implications 

CT. in RGWs detections and in cosmology. Analytical calculations show that, although 
the spectrum of RGWs on low frequencies is less affected by at ^ 0, but, on high 
frequencies, the spectrum is modified substantially. Investigations are made toward 

^ ■ potential detections of the at-modified RGWs for several kinds of current and planned 

^ I detectors. The Advanced LIGO will likely be able to detect RGWs with at > for 

2 ! inflationary models with the inflation index 13 = —1.956 and the tensor-scalar ratio 

O ' r = 0.55. The future LISA can detect RGWs for a much broader range of {at, j3, r), 

\ and will have a better chance to break a degeneracy between them. Constraints on 

^ at are estimated from several detections and cosmological observations. Among them. 



the most stringent one is from the bound of the Big Bang nucleosynthesis (BBN), and 



^ \ requires at < 0.008 rather conservatively for any reasonable {(3, r), preferring a nearly 

power-law spectrum of RGWs. In light of this result, one would expect the scalar running 
index as to be of the same magnitude as at, if both RGWs and scalar perturbations are 
generated by the same scalar inflation. 

PACS numbers: 04.30.-w, 04.80. Nn, 98.80.Cq 

1. Introduction 

Inflationary models predict a stochastic background of relic gravitational waves (RGWs) [1, 2, 3, 4, 5, 6], 
whose spectrum depends upon the initial condition when they were generated. After that, only the 
expansion of spacetime background will substantially affect its evolution behavior in a determined fashion, 
since the interaction of RGWs with other cosmic components is typically very weak. Therefore, RGWs 
carry a unique information of the early Universe, and serve as a probe into the Universe much earlier than 
the cosmic microwave background (CMB). 
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Unlike gravitational waves radiated by usual astrophysical process, RGWs exist everywhere and anytime, 
and, moreover, its spectrum spreads over a broad range of frequency, (10~^^ ~ 10^°) Hz. Therefore, RGWs 
are one of the major scientific goals of various GW detectors, including the ground-based interferometers, 
such as the ongoing LIGO [7], Advanced LIGO [8], VIRGO [9], GEO [10], and the space interferometers, 
such as the future LISA [11, 12], DECIGO [13], and ASTROD [14], the cryogenic resonant bar detectors, 
such as EXPLORER [15], NAUTILUS [16], the cavity detectors MAGO [17], the waveguide detector [18], 
and the proposed Gaussian maser beam detector at GHz[19]. Besides, a long period of observations of 
pulsar arrival times can be used as GW detectors at nanoHertz [20, 21], such as PPTA [22]. Furthermore, 
the very low frequency portion of RGWs also contribute to the CMB anisotropies and polarizations [23], 
yielding a magnetic type polarization of CMB as a distinguished signal of RGWs. WMAP [24, 25, 26, 
27, 28, 29, 30], Planck [31], and the proposed CMBpol [32] are of this type. Therefore, the detailed 
information of RGWs is much desired both for the detection of RGWs itself and for cosmology as well. 

The spectrum of RGWs depends on several physical factors. After being generated, RGWs will be 
affected by a sequence of stages of cosmic expansion, including the current acceleration [4], also by physical 
processes in the early Universe, such as the neutrino free-streaming [33, 34, 35], the QCD transition, and 
the e"'~e~ annihilation [36, 37], etc. But, over all, it depends most sensitively upon the initial condition, 
which includes the initial amplitude, the spectral index /?, and the running index at, as well. These 
parameters are predicted by specific models of inflation scenarios [4]. 

In this paper, based on our previous analytic work of RGWs [4, 35, 37], we will study the modifications 
by a non-zero at, together with upon the spectrum, and explore its implications for gravitational 
wave detections. By specifying the scale factor for various expanding stages and by a normalization in 
light of CMB temperature anisotropies of WMAP 5-year [28, 29, 30], the analytic spectrum h{k,TH) 
of RGWs with dependence on at and /3 will be demonstrated. Although a small at seems to have 
insignificant influences upon the spectrum of RGWs on very low frequencies, it will cause increasingly 
substantial modifications upon the spectrum on higher frequencies. This inevitably leads to far-reaching 
consequences to RGWs detection, since most of detectors are, or will be, operating at various medium 
and high frequencies, from around 10~^ Hz up to around 10^ Hz. The previously estimated constraints 
on RGWs should be revised in the presence of at accordingly. To this end, comparisons will be carried out 
between the theoretical spectrum of RGWs and the sensitivity of various ongoing and forthcoming GW 
detectors. Thereby, constraints on at will be derived and their implications in cosmology will be discussed. 

The outline of this paper is as follows. In section 2, the scale factor a(T) is specified for consecutive 
stages of cosmic expansion, and the construction is briefly reviewed for the analytical solution of the 
RGWs. In section 3, we present the resulting spectrum of RGWs with a scalar running index at and 
demonstrate the induced modifications. In section 4, comparisons are made between the calculated RGWs 
and the sensitivity of several kinds of ongoing and planned GW detectors, thereby, constraints on RGWs 
are obtained and implications in cosmology are discussed, in this paper we use unit with c = h = ks = 1- 



2. Analytical Solution of RGWs in Expanding Universe 



For a spatially flat {k = 0) universe the Robertson-Walker spacetime has a metric 

= a^{T)[-dT^ + dijdx'dx^], (1) 
where r is the conformal time, and the scale factor a(T) is determined by the Friedmann equation 

-P, (2) 
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where ' = d/dr. From the very early inflation up to the present accelerating expansion, a(T) can be 
described by the following successive stages [2, 4]: 
The inflationary stage: 

a(r) = Zo|r|'+^, -oo < r < n, (3) 

where the inflation index (3 is an important model parameter, related to the spectral index of primordial 

perturbation via = 2/3 + 5. The special case of /3 = —2 corresponds the exact de Sitter expansion. 
But both the model-predicted and the observed results, such as WMAP, indicate that the value of f3 can 
differ slightly from —2. In our presentation, (3 = —2.015 and P = —1.956, corresponding to Ug = 0.97 
and Ug = 1.089 respectively, are also taken for illustration. 
The reheating stage [35, 37]: 

a{T) = a^\T - Tpl^+f^' , Tl<T<Ts. (4) 

As a reheating model parameter, we will mostly take Ps = —0.3. 
The radiation-dominant stage : 

a{T) = ae{T - Te), Tg < T < T2. (5) 

The matter-dominant stage: 

a(r) = amir - Tmf , T2 < t < te- (6) 
The accelerating stage up to the present time th [4]: 

a{T) = Ih\t - Ta\'~'^ , TE<T<TH, (7) 

where 7 is a r^A-dependent parameter. For instance, 7 ~ 1.06 for Oa = 0.65, and 7 ~ 1.044 for 
i^A = 0.75 [35]. To be specific, we take 7 ~ 1.044 and J^a = 0.75 in this paper. 

In Eqs. (3) - (7), the five instances of time, ri, r^,, T2, te, and th, separate the different stages, 
and can be determined by the relations [35]: = 300 for the reheating stage, = W"^^ for the 

radiation stage, ^ = ^(pdntEl = 3454^^ for the matter stage, and ^ = for the 

=> a(T2) a(T2) a(TH) a(TH) * a(TE) ^ilm' 

present accelerating stage, and th is to be fixed by the normalization 

\TH-Ta\ = l. (8) 

In the expressions of a(r), there are twelve parameters, among which /3, (3s and 7 are imposed as the 
model parameters. By the continuity of a(r) and of a(r)' at the four instances ri, Tg, T2 and te, one can 
fix other eight parameters. The remaining Ih can be fixed by 



where Hq is the present Hubble constant. We will take the Hubble parameter — 0.71. Thus a(r) is 
completely fixed [35]. 

In the present universe the physical frequency for a conformal wavenumber k is given by 

The comoving wavenumber kn corresponding to a wavelength of Hubble radius I/Hq at present is given 

t, = ?^ = 2.-,, (11) 
and another wavenumber which will be used is 

kE ^ = kHin^/n^)^/'\ (12) 

i/H{TE) 

whose corresponding wavelength at the time te is equal to the Hubble radius 1/H{te) at that moment. 
Note that, in Eq. (12) we have made corrections to that in Ref. [35]. 
In the presence of the gravitational waves, the perturbed metric is 

ds'^ = a'^{T)[-dT^ + {Sij + hij)dx'dx^], (13) 

where the tensorial perturbation hij is traceless /?/• = and transverse hijj = 0. It can be decomposed 
into the Fourier /c-modes and into the polarization states, denoted by a, as 

h^Jir,^)= E j 4^'^^^k\ry''-'^ (14) 
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where h^^l*{T) = h^j^^r) ensuring that hij be real, efj is the polarization tensor. In terms of the mode 
h^]^\ the wave equation is 

h^:^"{r) + ^^h^:^'ir) + k^ht\r) = 0. (15) 

Assuming each polarization, x, +, /i^^^(r) has the same statistical properties, the super index (cr) can be 
dropped As listed in Eq.(3) through Eq.(7), the scale factor has a power-law form 

a(r) cx t"", (16) 

and the solution to Eq.(15) is a linear combination of Bessel and Neumann functions 

/ifc(r) = r^"[CiJ^_i(A;r) + C2iV^_i(A;r)], (17) 

where the constants Ci and C2 for each stage are determined by the continuity of and of /i^ at the 
joining points ti,Ts,T2 and te [4, 35, 37]. Therefore, the analytical solution of RGWs is completely fixed, 
once the initial condition during the inflation is given. 
For the inflationary stage, one has 

/ifc(r) = Aolo^lrl-^'^+'^^AiJi^^ikT) + ^2 J_(i+;3)(fcr)], -00 < r < n, (18) 



where the /c-independent constant Aq determines the initial amplitude, and 

Ai = l—.[le'-f^/\ A2 = iA^e-^-(', (19) 

cos /JTT V 2 

are taken [38, 4, 35], so that in the high frequency limit linijt^oo ^fe (''") e~'^^'^ the adiabatic vacuum is 
achieved [39]. In the long wave-length limit, kr ^ 1, the fc-dependence of /ife(r) is given by 

/ife(r) cx Ji+pikr) a A; 5+^. (20) 
3. Spectrum of RGWs with a Running Index 

The spectrum of RGWs h{k,T) at a time r is defined by the following equation: 

h\k,r)- = {0\h'^{x,T)hij{x,r)\0), (21) 



/o k 

where the right-hand side is the expectation value of the h^^hij. Calculation yields the spectrum as follows 

h{k,T) = ^k^/'\h{T)\. (22) 
TT 

Note that this expression has a factor \/2 in place the factor 2 in Ref. [35]. As the initial condition, the 
primordial spectrum of RGWs at the time of the horizon-crossing during the inflation is usually taken 
to be a power-law form [2, 4, 35]: 

hik,r,)=A{^f+^, (23) 
kn 

where the index /3 ~ — 2 for a nearly scale-invariant spectrum, and A is proportional to Aq in Eq.(18). 
In principle, both (3 and A are determined by the specific inflationary model. Here we take them as two 
independent parameters. In literature, the following notation is often used for the RGWs spectrum [24, 28] 

Al{k)^h\k,n) = Al{ko)I^^J , (24) 

where ko is a conformal pivot wavenumber, whose corresponding physical wavenumber is A:q = ko/a{TH)- 
For WMAP, the pivot k^ = 0.002 Mpc"^ is taken [24, 25, 28]. Comparing Eqs. (23) and (24) yields 

nt = 2/3 + 4 (25) 

and 

A = A„{ko) i^-^j . (26) 

For each cosmological model with kn and /? being given, Aft(fco) determines A. More often in literature, 
a tensor-to-scalar ratio r is introduced [24] 

where Ajj(A;o) is the amplitude of the curvature spectrum at k = k^, and has been fixed Ajj(A;o) = 
(2.41±0.11) X 10-9 by WMAP5 Mean [30], and A|(A;o) = (2.445±0.096) x 10"^ by WMAP5+BA0+SN 
Mean [28]. Note that this scalar amplitude has been "fixed" only after the assumption that the data do 



not contain gravitational waves, i.e. r = 0. Here we use r as a convenient representation of the amplitude 
normalization of Ah{ko) at ko, i.e., Ah{ko) = 4.94 x 10~^r^/^. Thus, Eq.(26) is rewritten as the following 

A = 4.94 X 10-^3 [hlLf+P, (28) 

At present, only observational constraints on r have been given by WMAP [28, 29, 30]. To be specific in 
our presentation, r = 0.55 and r = 0.22 will be taken, respectively. 

In general, the spectra of primordial perturbations, both scalar and tensorial, deviate from the exact 
power-law form except when the inflation potential is an exponential. As an extension, one usually consider 
the following form of power spectra [40, 41] 

/ U \ -H-ns(feo)+|as ln(fc/feo) 

A|(fc) = A|(fco) j , (29) 

Al{k) = Aliko) J-) , (30) 



which contain the "running" spectral indices ag = dug/dhik for the scalar perturbations and at = 
drit/dlnk for the tensorial perturbations. Currently, WMAP has given some preliminary constraint on the 
scalar index rig and the scalar running index a^. At the pivot wavenumber /cq = ko/a{TH) = 0.002 Mpc"^, 
WMAPl has given Us = 1.20l{]:J? and = -0.077t[J:J]^^ [24]; WMAP3 has given = 0.95ll[|:|ji^ and 
as = -0.055lo:o35 [26]. With a better determination of the third acoustic peak, WMAP5 has given 
an improved result: Ug = 1.0871^273 and = -0.050±0.034 [28]; and WMAP5+BA0+SN has given 
Us = 1.0891^;^^^ and = -0.05312:^28' °^ = 0.970 ±0.015 without [28]. Compared with n^, the 
value of the scalar running ag is relatively small. Thus WMAP5 data do not significantly prefer a scalar 
running index [30]. See also Refs. [42] for relevant discussions. 

But so far there is no direct observation of the tensorial index nt nor the running index at- In the slow 
roll inflationary models driven by a single scalar field, the tensorial indices, nt and at, are determined by 
the inflationary potential and its derivatives, so are the scalar ones, Ug and as, as well [40, 41]. There 
would be relations between the tensorial indices and the scalar ones, if one imposes further a consistency 
relation. In our context, for generality, we will treat nt and at as parameters independent of and ag. 
Corresponding to Eq.(30), the primordial spectrum in Eq. (23) is modified to 

h{k,n) = A{^f+^Aa,ik) (31) 
kn 

where the the extra factor 

^„,(A:) = (-^)>''^(W (32) 
ko 

is the at— induced deviation from the simple power-law spectrum, reflecting an extra bending. With the 
help of Eq.(28), this is 

h{k,Ti) = AjR(A;o)r^(A)2+/3^ (fc). (33) 

For a tiny at = 0.01, in the very low frequency range, from v = uq = 3 x 10"^^ Hz to i/ = 3 x 10~^^ 
Hz, Aat{k) goes from 1 to 1.055, causing only a minor increase in the amplitude by < 5.5%. However, at 
very high frequency, say u = 10^ Hz, Aat{k) ~ 10^, enhancing the amplitude by 4 orders of magnitude. 
Thi<; rlra<;tir pffprt rpniiirp<; a rlptailprl inup<;ticra1-inn inlri rvj. and il-<; rnn<;pniipn1-ial imnlira1-inn<; 




Figure 1: The spectral energy density ^gii') with = for two sets of r and /3. 



As discussed in Refs. [4, 35, 37], at the present time r = th, the long wavelength modes /iA;(r) with 
wavenumber k < kE are still outside the horizon, their spectrum are still of the form in (33), h{k,TH) = 
h{k, Ti). Given this initial condition with a running index, the analytic calculation of the spectrum of RGWs 
can be carried out straightforwardly, in the same way as the non running case [4, 35]. The only difference in 
the actual computing procedure is the amplitude normalization of the RGWs spectrum, which can be taken 
at the wavenumber kE, corresponding to a physical frequency ue = kE/27Ta{TH) = Hq/{1 + ze) ~ 10~^^ 
Hz. With the help of Eq.(33), it is given by 



h{kE,rH) = A^(A:o)ri(^)2+/^+i-'-fe/'=o). 

ko 



(34) 



In a cosmological model with a given set of {r,/3,at), the resulting RGWs spectrum h{v,TH) at present 
is fully determined. 

Another important quantity often used in constraining RGWs is its present energy density parameter 
defined by = where pg = 32^^ ^^i.o^^o energy density of RGWs, and = 3H§/8ttG is the 

critical energy density. A direct calculation yields [2] 



n 



gw 



ng{iy) — , 



with 



(35) 



(36) 



being the dimensionless spectral energy density. From this expression, it is seen that the spectral energy 
density 0,g{i>) and the spectrum h{v,TH) are two equivalent quantities, and h{v,TH) / is just the 
characteristic amplitude, denoted by hc{f) in Ref.[5]. As the cutoffs of frequencies, the lower and upper 
limit of integration in Eq.(35) can be taken to be uiow — 2 x 10^^*^ Hz and Vupper — 10^*^ Hz, respectively 
[35]. In the case of at = 0, the slope of 0,g{u) is fixed by the index (3, and the overall amplitude of 
Qg{v) is fixed by r, as shown in Fig. 1. In the following, two combinations (r = 0.55, /3 = —1.956) and 
(r = 0.22, /3 = —2.015) will often be taken for specific illustrations. 

When there is a running index at 7^ 0, the slopes of 0,g{v), and of h{u,TH) as well, are affected 

t;iih<;tan1-iallv Ficr 9 rlpmnn<;1-ra1-p<; h(i/ tzt^ fnr wariniit; usiiipi; nf n/j. in thp mnrlpl \A/ith r = f) 99 anH 




Figure 2: The spectrum of RGWs for various at- 

(5 = —2.015. It is seen that a greater at yields a higher h{v,TH) and that the modifications due to at 
increase with the frequency v. While the effects are small in the low frequency range, they are quite 
substantial in the high-frequency range. For instance, in going from at = —0.01 to at = 0.01, the 
amplitude of h{v, th) gets enhanced by 3 orders of magnitudes at z> ~ 10^^ Hz falling the range for LISA, 
5 orders at ~ 10^ Hz for LIGO, 6 orders at ~ lO'' Hz for MAGO [17], and 9 orders at ~ lO^Hz for 
the Gauss beam [19]. Equivalently, Fig. 3 gives the a^- dependence of the spectral energy density Q.g{v), 
which shows more drastically the variations in high frequencies due to at. 

Notice that the ratio r, the index (3, and the running index at in Eq.(31) play different roles in shaping 
the spectrum of RGWs. r sets the amplitude, [3 fixes the overall slope of the spectrum, and at gives 
an extra bending to it. However, in a rather narrow interval of detecting frequencies, r, /3, and at have 
a degeneracy to certain extent, since a larger value of each of them tends to enhance the amplitude of 
h{v,TH) in the interval. The narrower the interval is, the stronger the degeneracy will be. Therefore, if 
a detector is operating in a narrow interval of frequencies, it can only detect RGWs for a combination of 
parameters (r, /3, at)- Those detectors operating over a broad frequency interval will have a better chance 
to break the degeneracy. 

4. Constraints from Detections and Implications 

The a^-induced modification of the RGWs spectrum has practical implications for the ongoing and 
planned GW detections in the medium, and the high-frequency ranges. Some previously estimated con- 
straints on RGWs were based on the theoretical spectrum without a running index [4, 35, 37]. Now these 
will be subsequently revised to certain extent. 

Figure 4 gives the comparison of the sensitivity curves of LIGO [7] and Advanced LIGO [8] with the 
theoretical spectra of RGWs for various parameters (r, j3, at), in the frequency range (10^, 10^) Hz. Note 
that, in order to compare to the strain sensitivity hf{v) [5, 43] of the detectors, the amplitude per root 
Hz, h(y^Tii) I ^ , has been used. It is seen that, for r = 0.55 and /3 = —1.956, the LIGO I SRD [7] has 
already put a constraint on the running index: at < 0.013. It will yet not be able to detect the signals of 
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Figure 4: Comparisons of the spectra with the sensitivity of the LIGO I SRD Goal achieved by LIGO 
S5 [7] and of the Advanced LIGO [8]. 

RGWs for r = 0.55, j3 < —1.956 and at < 0.013. On the other hand, with a substantial improvement in 
sensitivity, Advanced LIGO will be able to detect the RGWs from models with r = 0.55 and /? > —1.956 
and at > 0, but still it will unlikely be able to detect RGWs for r = 0.22 and f3 = -2.015 even for 
at = 0.01 and less. 

Figure 5 is a comparison of the theoretical spectrum with the LISA sensitivity curve [11, 12] for the 
ratio S/N = 1 in the frequency range (10~^, 10*^) Hz, where one year observation time has been assumed, 
which corresponds to a frequency bin Au ~ 3 x lO^^Hz around each frequency. To make a comparison 
with the sensitity curve, one needs to rescale the theoretical spectrum h{u,TH) into the rms spectrum 

h{v,TH , IS.u) in the band Av [2, 5], 

h{u,TH,Au) = h{v,TH)'J—. (37) 

The plot shows that LISA by its present design will be quite effective in detecting the RGWs around a 
rancrp nf frpniipnrip<; (^(^~^ 15 x 10"^^ H7 In narl-iniiar the nint tplk thai- I ISA will he ahle fn dptert 
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Figure 5: Comparison of the rms spectrum with the LISA sensitivity curve [12]. 




Log„v [HZ] 

Figure 6: Comparison of the theoretical spectrum /i(z^, r//) with the converted sensitivity curve of 
LIGO, Advanced LIGO, LISA, and DECIGO [13]. The parameters of the four top spectra h{v^TH) 
are the same as the previous Fig. 5. 

RGWs for parameters r > 0.22, (3 > -2.015, and at > 0. Thus, regarding to detection of RGWs, 
LISA is expected to perform much better than LIGO and Advance LIGO. This advantage is due to the 
property that the RGWs have a higher amplitude in the range of lower frequencies. The situation is 
illustrated in Fig. 6, in which the sensitivity curves of LIGO and LISA are converted pertinently, in order 
to compare with the theoretical spectrum h{v,TH)- Moreover, regarding to the RGWs detection, the 
frequency range covered by LISA is very broad, as compared with LIGO. This is because the low frequency 
portion of LISA sensitivity curve has a slope that is rather close to that of RGWs spectrum in the involved 
region. This feature of LISA is important and can be instrumental in breaking the (r, /?, a^) degeneracy, 
as mentioned earlier. Besides, in Fig 6 the sensitivity of planned DECIGO [13] is also presented, together 
with two more spectra calculated for very low running indices at = —0.01 and at = —0.03, in the model 
(r = 0.22, /3 = —2.015), respectively. DECIGO, if implemented, will be a powerful detector, capable of 
detecting RGWs with a very low running index at > —0.03. 
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Figure 7: Comparisons of the spectrum with the EXPLORER sensitivity curve [15]. 
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Figure 8: Comparisons of the spectrum with the MAGO sensitivity curve [17]. 

Fig. 7 compares the theoretical spectrum hiv^Tu) I with the 2005 run sensitivity curve of cryogenic 
resonant bar detector, EXPLORER in the frequency range (890,920) Hz [15]. It is seen that, even for 
13 = -1.956 and at = -0.01, RGWs are still far beyond the reach of EXPLORER by three orders of 
magnitude. 

Fig. 8 compares the theoretical spectrum h{v,TH) / y/v with the sensitivity curve of MAGO, a double 
spherical cavity detector, around the frequency 4002 Hz [17]. For /5 = —1.956 and at = 0.01, RGWs are 
also beyond the reach of MAGO by three orders of magnitude. 

The proposed Gaussian maser beam detector will operate at a very high frequency, say ~ 4.5 GHz [19]. 
Still the spectrum of RGWs stretches to such high frequencies with enough power, depending on both 
the energy scale involved in the specific inflationary models and the reheating process. For the designing 
parameters for the proposed detector, its sensitivity is short by about 6 orders of magnitude. As for the 
prototype loop waveguide detector [18] operating at a high frequency 100 MHz, the sensitivity at present 
is short by much more. 

For lack of direct detection of RGWS, cosmological considerations can be more effective in providing 



the constraints upon the RGWs through its energy density. In particular, the BBN process occurring at 
a temperature T ~ a few MeV during the early universe is sensitive to the total cosmic energy density, 
including that of RGWs. An increase in the energy density of RGWs will enhance the freezing temperature 
of neutrons and hence the light-element abundances. Consequently, this will bring observational constraints 
on the running index at- The energy density of RGWs should not be too large, otherwise, it will significantly 
affect the outcome of the BBN process. Measured abundances of light-element constrain the number of 
additional relativistic species at BBN to an equivalence of SN^, = 1.4 neutrino degrees of freedom, 
corresponding to hQ^gyj < 7.8 x 10~^ [44, 5]. Besides, if the initial perturbation amplitude of RGWs 
is non-adiabatic, as is the case of RGWs being generated during inflation, measurements of the CMB 
power spectrum also provide a constraint hQQg{u) < 8.4 x 10~^ in the very low frequency range v = 
(10~^^, 10~^°) Hz [45], comparable to that from BBN. In constraining the energy density of RGWs, 
sometimes 0.g{ij) and $7^^^, were used interchangeably in literature. But it is at most an approximation, 
valid only under the condition that the integration interval du/u = dlogu ~ 1 and ^g{u) is rather flat 
and smooth. In this paper we distinguish and and give a more accurate treatment. Using 

the normalization of Eq.(34) and carrying out the integration in Eq. (35), we obtain the a^— dependence 
of the the energy density parameter Qg^ in Fig. 9 for various values of parameters. Adopting the BBN 
bound, the constraint on at is found to be at < 0.0014 for r = 0.55 and P = —1.956, at < 0.0056 for 
r = 0.22 and /? = -2.015, and at < 0.0077 for r = 0.001 and /? = -2.0, respectively. One can infer 
that the constraint at < 0.008 for any reasonable set of cosmological parameters. It should be mentioned 
that, although the BBN bound yields a constraint on RGWs, it does not provides a direct detection of 
RGWs, so one can only get upper bounds of at- 

This rather stringent constraint on at supports a nearly power-law spectrum of RGWs, and is consis- 
tent with the expectation from scalar inflationary models [40]. However, a comparison shows that, the 
magnitude of at constrained by the BBN bound is even smaller by one order than the scalar running 
index ag obtained by WMAP on large scales [24, 26, 28, 30]. If both RGWs and scalar perturbations are 
generated by the same inflation, one expects to be nearly as small as at for several kinds of smooth 
scalar potential [40]. In light of the stringent constraint on at, it is hinted that ag should also be rather 
small. There has been some debate on the significantly non-vanishing scalar running ag [42]. Therefore, 
it is much desired that constraints on the scalar running index be drawn from observations on smaller 
sales. 

Another kind of stringent constraints on RGWs come from observations of millisecond pulsars, which 
can serve as a gravitational wave detector in the low frequency range. For instance, by analyzing the 
uncertainty e in the arrival timing of pulses for a duration T of observation, the pulsar will be sensitive 
to gravitational waves h{u,TH) ~ e/T with u ~ 1/T. For PSR B1855+09, a bound has been given for 
RGWs [5, 21]: J7g(z/*)/ig < 4.8 x 10'^ {i^ / 1^^)^ for ly > z/*, where i^^ = 4.4 x 10"^ Hz. Another treatment 
gives r2g(i/*)/iQ < 2 X 10^^ {u / u^,)'^ for u > i/*, where z/* = 1.9 x 10^^ Hz [46]. Applying this bound to 
compare with the calculated energy density spectrum ^giiy), we obtain the constraint on RGWs, shown 
in Fig.lO. For the parameters r = 0.55 and (3 = —1.956, the pulsar detector puts a constraint at < 0.01, 
which is less stringent than the BBN constraint. As an extension of this technique, Parkes Pulsar Timing 




Figure 9: The at-dependence of energy density parameter for various r and (3. 




Figure 10: Comparison with the constraint from data of mihisecond pulsars [46] and PPTA [47]. 

Array (PPTA) [22, 47] consists of a sample of 20 millisecond pulsars distributed over the entire sky, 
and correlations in the time residuals of pulsars help to disentangle RGWs signals. After five years of 
observation, it will improve the sensitivity by about one order over that from a single pulsar and will have 
a chance to detect RGWs of r = 0.55, /3 = -1.956, and a* = 0. But for r = 0.22 and (3 = -2.015, the 
detector will still not be able to detect RGWs for a small value of at. 

In summary, by allowing for a tensorial running index at, the RGWs spectrum h{v,TH) will be sig- 
nificantly affected, particularly in higher frequencies. A positive at tends to bend the high frequency 
portion of h(y, th) to upward, while a negative at will do the opposite. This has brought about significant 
consequences to various ongoing and planned detectors, and led to reexaminations of the constraints on 
parameters of RGWs that were previously estimated for the case at = 0. For instance, a small variation 
of at from —0.01 to 0.01 will increase the amplitude RGWs by several orders of magnitude, depending 
on frequencies. It is interesting to note that LISA by its design will be able to detect RGWs with at > 
for the parameters r = 0.22, j3 = —2.015. For an inflationary model with r = 0.55 and j3 = —1.956, the 
observational constraint is at < 0.013 from LIGO S5, at < 0.01 from millisecond pulsar PSR 61855+09. 



The most stringent constraint coming from the BBN bound is at < 0.008 as a rather conservative esti- 
mate for any reasonable set of cosmological parameter r). The resulting tiny at prefers the simple 
inflationary models with a nearly power-law spectrum of RGWs, and would also hint a rather small scalar 
running index within scalar inflationary models. It is also found that there is a degeneracy of at with 
/? and r in a narrow interval of detection frequency. Detectors working in a broad frequency range, such 
as LISA, may have a better chance in breaking the degeneracy. 
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